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1. Summary of theoretical background 


This brief article gives an update on [j^ in the sense that we here analyze a very much enlarged 
new data set. This is found in the next section while we here summarize the theory for the reader’s 
convenience*. 

We consider the single component Z(2) symmetric scalar field theory on a torus of length 
L embedded in four dimensional Euclidean space. We employ the simplest hypercubic lattice 
discretization with L/a sites in each direction and the standard nearest neighbor lattice action 




-2ff£(p(x)(p(x + /t). 

X^i 


( 1 . 1 ) 


The standard picture for this quantum field fheory is |^, Q] fhaf fhere is a critical line fQ.(A) 
where fhe model possesses a continuum limif. This limif in fc may be approached from below 
fo reach fhe symmefric massive continuum fheory or (for A > 0) from above fo define fhe broken- 
symmefry massive model. The laffer is of particular fheorecfical inferesf here due fo some similaiify 
fo fhe Higgs field in fhe sfandard model. 

As a sfandard way fo renormalize fhe infinite volume fheory in fhe broken phase we may mafch 
fhe Fourier fransform G{p) of fhe fwo poinf correlafion funcfion 


X 


( 1 . 2 ) 


fo fhe asymptotic form 


G{p)=Z{\^d^{p) + 


p^ +rn^ + 0{p^) 


(p^O) 


(1.3) 


in fhe limif of vanishing p. Here fhe laffice momenfa implied by our discretization 

2 

p^ = - sin {ap^jl) 


(1.4) 


have enfered. In fhis formula v and m are fhe renormalized vacuum expecfafion value and fhe renor¬ 
malized mass and Z is a mulfiplicafive renormalizafion factor. We prefer fo define v from a zero 
momenfum confribufion fo fhe unsubfracfed fwo poinf funcfion rafher fhan from fhe direcf expec¬ 
fafion value {(p) because we fhus avoid sublefies wifh an ofherwise necessary symmefry breaking 
external field and we prepare fhe exfension of fhe scheme fo a finite size system. 

On a finife torus - laffice or confinuum - all momenfum componenfs gef quanfized fo integer 
multiples of 27r/L. We now focus on fhree momenfa wifh fhe smallesf mufually differing values of 


P = (0,0,0,0) 

(1.5) 

2ti 

p* = (1,0,0,0) — 

(1.6) 

.271 


p„ = (1,1,0,0) — , 

(1.7) 


* A more detailed elementary introduction can be found in [^. 
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for which we enforce as exact equality and then solve for Z,m, v. The result is 


= {niLf 


G jp**) pl^L^ -G{p*) pIL^ 
Gip,)-G{p,,) 


( 1 . 8 ) 


and 


= {vL) = 


G(0) 


1 


■ -z 


(1.9) 


G(p*) L'^pI + z^ 

where we have introduced the dimensionless finite size scaling quantities z = mL and w = vL. In 
addition we have substitued 


5\P) 




pfl 


( 1 . 10 ) 


for the finite size system. A renormalized coupling in the broken phase is conveniently defined by 
fhe rafio of mass fo expecfafion value, 


g = 




( 1 . 11 ) 


In an expansion around one of fhe degenerafe minima g is seen fo coincide wifh fhe bare coupling 
up fo loop correcfions. 

In our numerical invesfigafion we have resfricfed ourselves fo fhe limif A —)■ oo. Then fhe pafh 
infegral over laffice fields (p wifh weighf exp (—S') reduces fo fhe Ising model where we sum over 
cofiguralions {(p(x) = ±1}. 

The required observable G{p) can be very convenienfly esfimafed in fhe loop represenfafion 
[^] of fhe Ising model which is efficienlly sampled by fhe worm algorifhm [^]. In fhis ensemble fhe 
lanh(2ff) expansion of 

^ (u) (p (v) (1.12) 

w,v (p 

is sampled. As a consequence fhe disfribufion of u and v is related fo fhe fwo poinf correlation, 


{(p{x)(p{0)) 


{{Sxu^v)) 

((5„,v)) 


(1.13) 


where double angles refer fo fhe average defined by ( |1.12D . Finally fhe desired Fourier fransforms 
are given by 


Gip) = 


((e-'>(““'’))) 

((5„,v)) 


n^cos i^Pf, (m-v)^ 


(1.14) 


where fhe reflecfion invariance in each direcfion is used fo gef a real represenfafion in ferms of 
cosines only. 


2. New data 

In comparison fo [||] we have subsfanlially exfended our simulafions. In fable |i| we compile 
our complefe dafasef. Aparf from some memory opfimizafions, our implemenfafion of fhe worm 
algorifhm is a sfandard one. The only 4D field fhaf we keep in memory, is fhe link field fhaf 
represenfs a graph of fhe lanh(2fc) expansion. Since if can only assume fwo values per link, V /2 
byfes suffice for ifs sforage. Thus even our largesf laffices fif comforfably info fhe memory of a 
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Lja 

2 k 

z2 

g 

'?lz^=10 

8 

0.1524600 

10.024(96) 

29.13(30) 

29.70(26) 

12 

0.1509920 

10.008(98) 

24.88(26) 

25.09(22) 

16 

0.1504500 

9.928(80) 

22.30(19) 

22.44(16) 

24 

0.1500460 

9.974(98) 

19.65(21) 

19.70(18) 

32 

0.1498990 

9.970(65) 

17.93(13) 

17.96(11) 

48 

0.1497900 

10.189(53) 

16.121(93) 

16.074(77) 

64 

0.1497484 

10.03(11) 

15.03(18) 

15.03(15) 

80 

0.1497294 

10.10(14) 

14.36(22) 

14.34(19) 

160 

0.1497035 

10.30(38) 

12.06(49) 

12.01(42) 


Table 1: Simulation results for z and g at growing Lja and for g corrected to refer to = 10. 


standard desktop PC. It took roughly 43k core hours to generate our most expensive (Lja = 160) 
ensemble. This corresponds to 1.2 x 10’ X W worm updates. 

The coupling g is related to the renormalized coupling g by 

§ ^T\ 2 2 /1 fXl ’ 

i. e. it only differs by a small lattice artefact. The rationale 0 is that the Callan Symanzik p 
function for g has a tree level artefact contribution which is absent for g. For our z’ = 10 data this 
amounts to a relative 1.25L^^ correction that is completely insignificant except for the smallest 
lattices. 



Figure 1: Data and linear fit for g ^dgjdz- 

We had to tune fc to approach z^ = 10. Table |T] shows that we were often successful within our 
small statistical errors. To implement the remaining tiny correction leading to the last column we 
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have numerically determined the derivative dg/d{z^). This is relatively easy by using the relation 

( 22 ) 

which holds in the loop representation |^. In this formula t = tanh(2fc) is used, A can be any 
fc-independent observable and K is the total number of links occupied by lines. We see in figure 
[I] that beyond L = 64 the connected correlation ( |2.2[ ) is too noisy to get a signal and we had to 
extrapolate with the shown linear fit. Its form is suggested by perturbation theory. We emphasize 
however that any en'or in this procedure only affects a systematic correction in the final results that 
itself is only of the order of the statistical error. 


L/a 



Figure 2: Cutoff dependence of the coupling g. 

Our main result is now represented by figure ^ The curves show the evolution with the per¬ 
turbative renormalization group at 1,2 and 3 loop^ precision where the evolution is started (in both 
directions) from our most precise data point at L/a = 48. The shaded bands represent the small 
error of the initial value. We see a perfect match of all our data points with the three loop evolution 
and thus complete consistency with the (logarithmic) triviality scenario. 

^We have to note that the three loop term is taken from and refers to z = oo. According to experiences in the 
symmetric phase |Q] the value for = 10 is expected to be very similar. 
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